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It is discussed how the equilibrium properties of the Ising model are described by an Hamiltonian 
with an antiferromagnetic low temperature behavior if only an heat bath dynamics, with the char- 
acteristics of a Probabilistic Cellular Automaton, is assumed to determine the temporal evolution 
of the system. 
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I. INTRODUCTION 

In this note we discuss the equilibrium properties of Probabilistic Cellular Automata (PCA) reversible with respect 
to a Gibbs measure derived by a suitable Hamiltonian. A PCA is a lattice model with discrete variables which 
are subject to a probabilistic simultaneous updating in discrete time steps: all configurations are accessible in a single 
updating. 

PCA arise as an extended definition of Deterministic Cellular Automata in which the updating follows a set of 
deterministic local rules. The huge number of possible deterministic (or probabilistic) rules makes the topic of 
Cellular Automata overwhelmingly abundant. One of the most famous Cellular Automata systems is the Conway's 
"Game of life" Q]; in spite of the very simple deterministic majority rule assigned, the system, which is a kind of spin 
lattice, presents an extremely rich and complex evolution pattern. 

PCA have been studied in a wide variety of contexts, ranging from biology to the theory of automation, while their 
role in statistical mechanics has been thoroughly investigated only relatively recently in 1^-^] . The particular family 
of automata we study is obtained by implementing in parallel fashion the heat bath dynamics In other terms, we 
define a rule for the transition probabilities such that all single spins of a lattice arc updated simultaneously with heat 
bath rates. This amounts to define a Markov chain for the evolution of the spin system, having the characteristics of 
a PCA. 

We observe that the way of implementing the heat bath dynamics reflects into a qualitative modification of the 
equilibrium properties of the model. In particular, an Ising-like ferromagnetic Hamiltonian with two body interactions, 
defines a PCA reversible with respect to a Gibbs measure determined by an Hamiltonian allowing a low temperature 
antiferromagnetic phase. This behavior is absent if a serial dynamics is implemented, for which at most one spin of 
the system is updated at any time. 

The paper is organized as follows. In Section [n| we define the PCA under consideration and the spin model 



Hamiltonian, H, determining the heat bath single spin rates. In Section III the structure and the low temperature 



antiferromagnetic properties of H' are discussed. There follows a discussion of its equilibrium properties and ground 



states. Finally, in Section IV, it is considered the case in which the heat bath spin rates are determined by a (two 
body) next-to-nearest neighbor interaction. The dynamical generation of antiferromagnetic couplings is reviewed in 
the conclusive remarks. 

II. COUPLING PROLIFERATION IN REVERSIBLE PCA'S 



Let A be a finite two-dimensional square lattice and |A| its cardinality. For each x = [xi,X2),y = (2/1,2/2) £ Awe 
denote by \x — y\ the Euclidean distance on the lattice. Let a(x) E {— 1, +1} a spin variable associated to the site 
ieA; the space {1, — 1} A of configurations is denoted by S. 
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Let us consider a generic Hamiltonian H : a g S — > if (cr) € R and the corresponding equilibrium Gibbs measure 

e -PH(*) 

with /? the inverse of the temperature. We now define the heat bath single spin rates: given the site x G A, we consider 
the Gibbs equilibrium measure for a x with respect to a fixed configuration cr on A \ {x}. Letting a G {— 1, +1}, we 
have 

/ I \ exp{-/3ff(q,cr)} 

Px{l> ~exp{-/3H(a,a)} + exp{-pH(-a,a)} ' lj 

where (±o, cr) are the configurations equal to cr on A \ {x} and to ±a on x. Note that the normalization condition 
Pxip\a) + p x (—a\a) = 1 is trivially satisfied. 

We can now implement the heat bath dynamics in a serial fashion, namely we can consider the Markov chain at (t 
being the discrete time temporal variable), with transition probabilities 

p/ s _ J (l/|A|)p x (r?(x)|cr) if 3x G A such that a = r\ on A \ {x} 

^ a > r D-\ otherwise (6) 

for all cr, i] G S. The transition probabilities (||) are reversible with respect to the Gibbs measure (Q), i.e., the detailed 
balance condition is satisfied or, equivalently, the equilibrium measure is the Gibbs measure (|l|). 

A different point of view can be taken ||: we define the transition probabilities P(cr, rf) in such a way that all the 
spins are simultaneously and independently updated, in a parallel fashion, with the heat bath rates (|^). Thus, instead 
of (||) we consider the Markov chain cr t defined by 

P(o-,v) = Y[px(v(x)W) Vcr, 7? G S . (4) 

xeA 

This amounts to define a Probabilistic Cellular Automata (PCA). 

In general the equilibrium properties of the Markov chain (|4|) are not trivial, for instance it is not obvious that 
there exists a Gibbs measure such that the detailed balance principle is satisfied. 

Let us consider, now, the case of the two body interactions and suppose that the Hamiltonian has the form: 

H ( CT ) = ~ J ^,y cr ( ;E ) cr (y) _ /l ^ cr ( x ) > ( 5 ) 

l,t/£A J)£A 

where J x _ y G R are the pair couplings between spins at sites x, y, and h x G R is the external magnetic field acting on 
the spin at site x. For physical reasons we suppose that the pair couplings are symmetric, namely J XiV = J y ^ x for all 
x, y G A. The heat bath single spin rates are given by 

Px(a\a) = — — — j , , l r rpr = -[l + atanh/^O)] , (6) 

1 + exp {— p [H (— a, a) — H(a, a)\\ I 

where 

s x(v) = ^2 Jx, v o-{y) + h x , (7) 

y£A\{x} 

for any a G S and x G A. It is easy to show Q that the Probabilistic Cellular Automaton (||) with single spin rates 
([|) is reversible with respect to the Gibbs measure [i! on S associated to the Hamiltonian 

H V) = ^ log COsh [/^(cr)] . (8) 

i£A iEA 

In other words the detailed balance condition 

P(cr,7 ? )exp{- J ff , (cr)}=P( J7! cr)exp{- J H"( ? 7)} , (9) 

is satisfied for any a,rj G S. This means that H'(o~) is the equilibrium Hamiltonian of a system governed by H(o~) 
and evolving with the law (Q). 
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The choice of the two body interaction in equation (|5|) is strictly connected to the reversibility of the resulting 
Probabilistic Cellular Automaton ([!]) @,S). For example, consider the three body interaction Hamiltonian H(a) = 
J2x y zeA Jx,y,x^{x)a{y)a{z) with the three body couplings J x ,y,z symmetric with respect to permutations of the 
indices and such that J x , y .z 7^ if and only if x 7^ y =/= z 7^ x. Then the problem of showing the reversibility of the 
parallel dynamics can be reduced to the problem of finding a function : a G S — ► 4>{<j) G K such that 

00) - 3/3?y(x) ^2 J x,y,zcr(y)cr(z) = <f>(rj) - 3(3a(x) ^ J x , VtZ -q{y)-q{z) , 

y,z£A\{x} y,zeA\{x} 

which has no solution. 

Let us now discuss the main feature of reversible heat bath derived probabilistic automata. As it has been seen 
above, if the starting Hamiltonian is given by (||), then the Markov chain is reversible with respect to the Gibbs 
measure with Hamiltonian H' given by (||). It is clear that new kind of interactions, different from the one present in 
the original Hamiltonian H, arise when H' is considered. 

Suppose, for instance, that the starting Hamiltonian has range r > 0, namely J x ^ y = for any x,y € A such that 
\x — y\ > r. Then we have S x (a) — J2 y - o<\x-y\<r Jx,y&(y) + h x , a sort of average of the spins inside a ball centered 
at site x with radius equal to r. Hence, by expanding H' as a sum of potentials we will get all the possible couplings 
inside the ball, starting from the two body up to the N(r) body interaction, with N(r) the number of sites inside the 
ball. In some sense these new couplings are dynamically generated. In the following we will discuss few interesting 
particular cases. 



III. NEAREST NEIGHBOR ISING MODEL 

Let us consider the standard nearest neighbors Ising model with no external magnetic field, namely we consider (|^) 
with J XiV = J/2 for any x, y £ A such that \x — y\ = 1, J XiV = otherwise, and h x = for any x € A. The Hamiltonian 
H 1 is the sum of averages performed over the four site crosses centered at each site of the lattice. We then expect all 
the possible interactions inside the cross. 

As it as been seen in jL0| it is possible to extract the potentials and rewrite the Hamiltonian H' in the following 
way 

H'(a) = -J, vfrMv) "^^ <r(x)a(y) - J 3 ^ a(x)a(y)a(w)a(z) , (10) 

where the three sums (see Fig. la) are respectively performed over the pairs of next to the nearest neighbors (sites 
at distance V%), the pairs of third neighbors (sites at distance 2), the four site diamond shaped clusters (plaquettes 
with side length equal to \/2). The coupling constants are given by 

J 1 = Jlogc OS h(2^)~i/3J ) J 2 = ij l5 ^3 = ^log^^~-^J , (11) 

where "~" means the limiting behavior for (3 — > 00. There exist several possible ways to extract the potentials. A 
very natural one, in the case of spin variable, is to consider the function ip x (a) = logcosh[(/3J/2) X^gA <J {v)\i where 
^x = {y G A : \y — x\ = 1} is the set of nearest neighbors of site x, and its expansion 

Vx {a) = ]T c(X) [] <r(y) , 

XCA X yGX 

with the coefficients c(X) given by 

<re{-i,+i} A =: 2/ex 

where, we recall, |A X | = 4 is the cardinality of A x (sec JTT|). 

It is important to remark that the second nearest neighbor interaction, Ji, is positive and dominating; hence we 
expect a low temperature antifcrromagnctic phase to exist. What appears very remarkable is that we have derived 
an antiferromagnetic behavior in a purely dynamical way as a result of the coupling proliferation. If a parallel heat 
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bath Ising dynamics is implemented, the equilibrium Gibbs measure shows a low temperature antiferromagnetic phase 
despite the simple physical ferromagnetic coupling of the Ising model. This phenomenon is. obviously, absent if the 
Ising heat bath dynamics is implemented in a serial fashion. 

The equilibrium properties of the model can be understood by remarking that two independent models are found 
if the lattice is partitioned into two square sublattices with step y/2 (the even and the odd sublattice). Each model 
is, indeed, an eight vertex model with nearest neighbors coupling Ji, next to the nearest neighbors coupling J? and 
plaquette interaction J3. This model has been widely studied both in two |l2}-|l4|] and three |l5],[l6| dimensions. From 
( |TT| ) and the very well known properties of the two-dimensional eight vertex model we have that on each sublattice 
there are two coexisting low temperature phases, respectively with positive and negative magnetization. Hence, by 
combining in all the possible ways the two phases we get, for the original model, the 3 different low temperature 
phases corresponding to the three ground states -0p, ipi,ip2 (see Fig. 2). 



It is of some interest a direct study of the hamiltonian (10): ground states can be defined as those configurations on 



which the Gibbs measure associated to the Hamiltonian H' , is concentrated when the limit (3 — > 00 is considered, 
namely as the minima of the energy 

E(a)= lim ^P- = -J2\S x (a)\ , (12) 

uniformly in a G S. It is rather clear that with periodic boundary conditions there exist three coexisting minima 
ipo, ipi, ip2 S S (see Fig. 2), with energy — 4|A|, such that ipo(x) = +l s ipi(%) = (— l) Xl+X2 and ip2{x) = —1 for all 
x = (xi,x 2 ) <G A. Notice that ipi is the chessboard configuration. 

The problem is, now, to understand if this coexistence of different states persists at a finite small temperature, 
namely if the system undergoes a low temperature phase transition. We give an heuristic argument: at finite tem- 
perature ground states are perturbed because small droplets of different phases show up. The idea is to calculate the 
energetic cost of a perturbation of one of the four coexisting states via the formation of a square droplet of a different 
phase. If it results that one of the three states ipo,ip\,ip2 is more easily perturbed, then we will conclude that this is 
the equilibrium phase at finite temperature. 

A simple calculation, sec 0], shows that the energy cost of a square droplet of side length n of one of the two 
homogeneous ground states plunged in one of the two chessboards (or vice versa) is equal to 8n. On the other hand 
if an homogeneous phase is perturbed as above by the other homogeneous phase, or one of the two chessboards is 
perturbed by the other one, then the energy cost is 16n. 

Hence, from the energetical point of view the most convenient excitations are those in which an homogeneous phase 
is perturbed by a chessboard or vice versa. Moreover, for each state f/'Oi ^Pii ^2 there exist two possible energetically 
convenient excitations: there is no entropic reason to prefer one of the four ground states to the others when a finite 
low temperature is considered. This remark strongly suggests that at small finite temperature the three ground states 
still coexist. 



IV. NEXT TO THE NEAREST NEIGHBOR ISING MODEL 



Let us consider the Ising model with no external magnetic field and next to the nearest neighbor interaction, namely 
we consider (|j) with J XAJ = J/2 for any 6 A such that \x — y\ — \/2, J XiV = otherwise, and h x = for any 
x e A. 



It is possible to extract the potentials as seen in Section III. The Hamiltonian H' can be written as 



H'(a) = -J, v&My) - ff (*MiO - J 3 £ a(x)a(y)a(w)a(z) , (13) 

where the three sums (see Fig. lb) are respectively performed over the pairs of third nearest neighbors (sites at 
distance 2), the pairs of sites at distance 2y / 2, the plaquettes with side length equal to 2. The coupling constants are 
still given by (|ll|). 

In order to study this model we remark that if the lattice is partitioned into four square sublattices with step 2, 
then we obtain four independent models one on each sublattice, each model being again an eight vertex model with 
nearest neighbors coupling J 1; next to the nearest neighbors coupling J 2 and plaquette interaction J 3 . Hence, on each 
sublattice we have the two degenerate ground states with all the spins, respectively, equal to one and minus one. By 
combining in all the possible ways these two states we get, for our model, 2 4 = 16 different ground states. On the 
torus, namely when periodic boundary conditions arc considered, some of the ground states are equivalent, so we get 
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the seven states -^V ipi,...,ipe in Fig. 2 and 3. The fact that the phase transition persists at finite small temperature 
is, as in Section [II a straightforward consequence of the known behavior of the eight vertex model. 

In this note we observed the relation between the two body ferromagnetic interaction of the Ising model and the low 
temperature antiferromagnetic behavior of the equilibrium Hamiltonian obtained evolving the initial system with a 
parallel heat bath dynamics. Any lattice system with two body interactions and having a self-organization resembling 
the parallel dynamics here described in the simple Ising model, could evolve towards unexpectedly complex equilibrium 
states. The influence of the parallel dynamics on the physical properties of the equilibrium measure emerges as an 
interesting new feature that deserves further explorations. 
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Figure captions 



Fig. 1: Couplings J\, Ji and J3 for the hamiltonian ( pL 0|) (resp. (13)) are shown in (a) (resp. (b)). 
Fig. 2: The three ground states ipo,ipi,ip2 depicted from the left to the right. 
Fig. 3: The four ground states V^., V'S) V'e depicted from the left to the right. 
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